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Abstract 

We study the diffusion equation on the surface of a 4D sphere and obtain a Kubo formula for 
the diffusion coefficient. 
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Since numerical simulations of fluids or plasmas within hyperspherical boundary condi- 
tions are appealing, it is not without interest to dispose of an explicit expression for the 
diffusion coefficient D in terms of time averages of microscopic quantities.— Our starting 
point will be the diffusion equation 

Pp(M, = - ^A^^) P(M, t) = , (1) 

where the point M lives on the surface of the four dimensional (4D) sphere 1S3 of center O 
and radius R (a hypersphere for short). D{> 0) is the diffusion coefficient and A*^^ is the 
Laplacian appropriate to the geometry.-''^ Eq. (P) was solved recently by Nissfolk et ai; their 
expression for p(M, t) is however akward and does not allow to obtain an explicit expression 
for D. A more convenient expression for p{M,t) is obtained in this paper which yields an 
explicit expression for D. 

We first note that the Brownian motion of a particle in S3 may also be viewed as the 
Brownian rotation of a linear AD molecule. As the expression of the rotational diffusion 
coefficient of a 3D rotator is known since the work of Debye^^ it should be easy to extend 
this old result to the AD case. Following Berne and Pecora^ we first expand p(M, t) on the 
complete basis set of hyperspherical harmonica^ 

00 

p(M,t) = 5^5^pz.,c.(t)lL,«(0, (2) 

L=0 a 

where the unit vector ^ = OM/R denotes the orientation of vector OM in and Yi^ai^) is 
a AD hyperspherical harmonics. We recall that the quantum number L is a positive integer 
and that a = {1711,1712) where mj {i = 1,2) takes the L + 1 values nii = —L/2,—L/2 + 
1, . . . , L/2. The AD harmonics YL^a{$) satisfy the following properties?^ 

• (i) orthogonality 

j dn{^) yLj^)yL',c.'{^) = h,L' , (3) 

where dVt{$^) denotes the infinitesimal volume element (or AD infinitesimal solid angle) 
on the hypersphere of unit radius. 

• (ii) completeness 

Y.Kc.mLAO = ^''{i,i'). (4) 



where S^^ is the Dirac distribution for the unit hypersphere defined as 

J dni^)fi^)s'^^ite) = fi^'). (5) 

• (iii) addition theorem 

T.'^um.Ae) = pl{^ ■ e) , (6) 

where the dot denotes the usual AD scalar product, i.e ^ ■ = cos{ip) where ip is the 
angle between the two unit vectors ^ and ^' and 

n / L + 1 sm(L + 

Picosv^ =^r^ ■ f,. ■ 7 

Note that Rip is the geodesic length between the points M and M' of 1S3. 

Moreover the Y^ ai^) are the eigenvectors of the Laplace-Beltrami operator A'^^ with eigen- 
values —L{L + 2)/R'^!^lt follows from this important property that eq. (Q) is equivalent to 
the set of equations 

d L(L + 2)\ , , 
the solution of which reads obviously as 

PLAi) = PlAO) exp {-DL {L + 2) t/R^) . (9) 

We shall denote by p(M, t|Mo, 0) the solution of (H)) corresponding to the initial condition 
p(M, 0) = 5'5^(M, Mo) = R~^ 5^-'{i, ^0), i-e. the solution of 

Dp(M,t|Mo,0) = 5{t) R-' 6^'{t^o) ■ (10) 

It follows readily from eqs. and (jH) that the Green function p(M, t|Mo, 0) can be expressed 
as 

p{M,t\Mo,0) = (t < 0) 

p(M,t|Mo,0) = -l^r2,«(6)^L,«(0 exp{-DLiL + 2)t/R') (t > 0) . (11) 

L.ct 

Eq. (fTT|) can be further simplified with the help of the addition theorem which yields our 
final result 

p{M, t\Mo, 0) = ^Y1 ^^(^0 ■ ^) (-^^ + 2) V^') {t>0). (12) 



Some comments are in order. 

(i) The solution (fT^ is invariant upon rotations about the axis as expected. Eq. (fT^ 
is a trivial generalization of Debye's resullj^^ (for the 3D rotor the Tchebychev polynomial 
Pl in eq. has to be replaced by a Legendre polynomial, up to a multiplicative constant) 

(ii) Note that 

j R^dn{^) p{M,t\Mo,0) = 1 it>0) , (13) 

i.e. the probability is conserved, the particle does not evaporate. 

(iii) We have also 

hm p(M,t|Mo,0) = -4^, (14) 

i.e. the solution of the diffusive process is uniform after infinite time (the volume of is 
equal precisely to 27i^R^) 

(iv) Let us define, following Berne and Pecora^, the time correlation function 

Cdt) ^ (PL(^(O)-I(t))) 

= / / ^'^^(^) ^^(^0 ■ ^) PiM,t\Mo,0) , (15) 

where we have averaged the initial position uniformly on 1S3. As a consequence of the 
orthogonality properties of the Pl we find that 
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CL{t)= ' 7/ exp(-DL(L + 2)t//?^) . (16) 



Defining now the reorientational time tl of our AD rotor as 

Cl{0) 



00 







Tl = dt 

(17) 



DL{L + 2) ' 

we have now at our disposal an explicit expression of the diffusion coefficient D. Note the 
aesthetic relation ti/tli = L'{L' + 2)/ L{L + 2) which generalizes Debye's relation to the 4D 
rotor. 

(v) It does not seem possible to obtain from eqs. (|TT)|) an expression for D in term of 
the time autocorelation function of the velocity of the particle. Moreover one can deduce 
from eq. (fT^ that the mean square displacement of the stereographic projection of point M 
vanishes for t ^ 00. 



It remains now to show that the expression (fT^ of p{M,t\Mo,0) is equivalent to that of 
refO This can be done as follows. Let us rewrite the reduced p = 27r^i?^p(M, t|Mo, 0) as 



where K = Dt/B?. A priori the angle is in the range (0,7r), however since the function 
is formally even in ip we define p(— V^) = p('?/') for negative angles. This gives us a periodic 
function of period 27r defined for all G M. We introduce now the periodic function 

F{ip) = / d^' sin(^') (19) 
Jo 

which can be rewritten after some algebra as 



where Fq is some unessential constant independent of angle ip. At this point we recall 
Poisson summation theorem which states that for any function ip{x) holomorphic in the 
strip —a < '^z < a one has 




L=0 



(18) 




(20) 




(21) 



where 




(22) 



Applying Poisson theorem for the Gaussian we get 




(23) 



which, after differentiation yields for p 



which coincides with the result of ref|2| apart the prefactor which is not specified. 




(24) 
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